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ABSTRACT 


This report deals with some of the existing 
methods to "globally" approximate the roots of 
transcendental equations namely, Graeffe's method? 
Summation of the reciprocated roots, Whittaker-Bernoulll 
method and the extension of Bernoulli’s method via 
Koenig's theorem. The Aitken's 8 a process is used 
to accelerate the convergence. Finally, the suitability 
of these methods is discussed in various cases. 
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NOMENCLATURE 


C.P*U* = a measure of computer time 

M = a constant representing the upper bound for a set 

- a function analytic at Z = r 

h = values derived recursively from Bernoulli’s difference equation 

y “ the maximum modulus of i>(z) 

H v p “ a sequence of Hankel determinants 

p = radius of convergence of a power series 

k = a normalizing polynomial 

vP 

K -a Hankel determinant of order p 

vP 

^ = the remainder for the Taylor series after v terms 

P ™ an approximate to the root of least modulus using the method 

of summing the reciprocated roots* 


v 



CHAPTER I 


INTRODUCTION AND SUMMARY 

In the area of applied mathematics, there is an urgent need to 
effect more efficient means of approximating the roots of transcendental 
equations. Existing standard methods employed in finding the roots of 
algebraic equations are of little aid in obtaining satisfactory results 
in this situation. Incidently, there exist no formulated method for 
finding the roots of algebraic equations above the fourth degree [14] . 
The same fact applies to a transcendental equation when it is repre- 
sented by a power series that has been truncated above the fourth degree 
term. 

There have been various attempts made in locating or approxi- 
mating the roots of certain types of transcendental equations. One of 
the more basic means of obtaining information about the nature of the 
roots of an equation (algebraic or transcendental), such as, h(x)= 
g(x)-f(x)=0, is through graphing the functions g(x) and f(x) on the same 
coordinate axes and noting the point (s) of intersection. After having 
made some suitable approximation to the desired root, then an iterative 
method may be used to refine the approximation. 

On the other hand, there are methods which are capable of 
yielding, in a more consistent manner, information about the roots 
of a given transcendental equation. One such method is the Graeffe 
method [151. Graeffe' s method guarantees convergence to a root 
through repeated root squaring [4]. 

There are other methods, though not discussed in this paper. 
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that are ’self starting’ or 'global' in the manner in which they 
approximate the roots to transcendental equations. These methods 
include Rutishauser ’s q-d algorithm and Bairstow’s method. The 
above methods are fairly effective in determining the roots of 
polynomial equations, but posses certain inherent disadvantages 
such as their sensitivity to truncation error where power series are 
involved. Furthermore, they are impractical from the standpoint of 
their compatability with the digital computer where C.P.U, time 
may be crucial. 

There exist another method of great prominence and it is 
known as Bernoulli's method. However, its convergence is somewhat 
less rapid than Graeffe’s method. In addition, it may be more 
difficult to deal with complex roots [4]. 

Later, Koenig proposed a theorem which extends Bernoulli's 
method to non-algebraic equations [5]. In this extension, provisions 
are made for treating equations having "equal" roots or in the case 
in which there are complex conjugates. Besides computing the root of 
least modulus within the region of analyticity, it is theoretically 
possible to compute all the roots of a transcendental equation that 
are within the radius of convergence of its power series representation^] . 
However, in some cases the convergence of Bernoulli's method may be 
somewhat less than ideal and in this case one may hope to employ some 
method of accelerating this convergence. One method of achieving this 
aim is the use of Aitken's § process [5], 
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Another method, although less popular, can be derived from 
the Bernoulli method is the Whittaker expansion [3,11]. This method 
relates the coefficients of a polynomial or a truncated power series 
to the roots via determinants. 

Among the miscellaneous methods, the root of least modulus 
for a transcendental equation may be evaluated by utilizing the 
standard relationship that exist between roots and coefficients of 
the truncated power series representing a particular transcendental 
equation. 

Thus in chapter 3 these methods will be applied to certain 
transcendental equations in order to obtain information (on a 
comparative basis) regarding their suitability to the computer, speed 
of convergence, ease of computations, and accuracy of approximation. 

We will study the nature of the non-algebraic equation by 
studying the coefficients of its power series representation. From a 
practical standpoint, we? cannot consider all the terms in the series, 
hence, we must truncate it at some point where the truncation error [15] 
will be minimal. Afterwards, we consider the truncated power series as 
a special case of a polynomial. However since the original equation is 
likely to consist of meromorphic functions, we must have some previous 
knowledge about the nature of their convergence in the region of analyticity. 
The following well known theorems and definitions are relevant to our 
further discussion in later chapters. 
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Roots and coefficients. Given the polynomial equation 

a 0 x n + a^x 11 ”^ + ... + a n = 0 , 

the following relationships exist between coefficients and roots [14]: 

(a) aj/a o - + r 2 + tj + ... + r n ) 

(b) a 2 /a 0 = (r 2 r 2 + r x r 3 +v.s 4- ri r n + ... + r n _ 1 r n ) 


(c) a 3 /a 0 = ~(f 1 r 2 r 3 + + . 


+ r ,r_ ,r_) 
ti“*z n-_L n 


(d) a /a - (-l) n (r r r ...r ) 
no 1 2 3 n 


Theorem 1. Every power series % c n (z-z o ) n has a "radius of convergence" 
R such that when 0 < R < the series converges absolutely for 
jz - z Q ;| < R and diverges for Jz - z Q | > R. When R = 0, the series 
converges only for z = z q . When R = the series converges for all z. 
The number R is given by 

R = 1/ lim yjc n | 


Theorem 2. Let f(z) be analytic in the interior of a circle C with 

center at z and radius R. Then at each point z interior to C 
o 

f(z) = Sf^(z )/n!(z-z ) n 
o o 

Definition 1. A function f(z) which is analytic in a domain D, except 
at some points of D at which it has poles , is said to be mememorphic. 



■CHAPTER II 


THEORETICAL DEVELOPMENTS OF GLOBAL METHODS 
Among the existing methods of approximating the roots of poly- 
nomial equations, Bernoulli’s method [1] is perhaps among the oldest kmown. 
In its original form, Bernoulli’s method was proven valid for algebraic 
equations* By expressing the equation in terms of a power series, it 
would seem possible that a modified Bernoulli method could be used 
to determine information concerning non- algebraic equations. 

Due mainly to the efforts of Koenig , the Bernoulli method was 
shown to be valid in the case of non-algebr.aic or transcendental 
equations. The theorems which follow were reproduced in order to 

(1) extend this method to the case of transcendental equations and 

2 

(2) show that Aitken’s 6 process can be applied jointly with the 
Bernoulli method. 

THEOREM A (Koenig’s theorem) Let h(z) = g(z)/f(z) - c Q + c-^z + ... be 
convergent for lzt<(R where f(z) has a simple zero at z « r and 
g(z) be any function which is analytic throughout the circle 
but g(r) i 0. If !r|< 5R<R, then c v /c y+1 = r + 0(a V+1 ) [5]. 

From the hypothesis, since z “ r is a pole then it is apparent that 
^(s) = (r-z)(h(z)) = b Q + b^z + b + ... is analytic for |zi<(R. 

Let ir| ( SR< R and v"'|r| < p < R, then ^(p) also converges having a term that 

a. -y 

is maximal in modulus which we call y . Then for any v , ib v |^yp and 
this establishes the fact that each is bounded. Expanding y (z) and 
equating similar coefficients, then 

* ’ ' : “ 

See [14] for a detailed development of Graeffe’s method. 
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rc.. - c , = b 

v v-1 v 

Upon multiplying the first v+ 1 equations by 1, r, r 2 , r^, . ,.,r v 
and adding we have r = ^(r) where^ v (z) = b Q + bj^z + ... + b v zV 
Letting R y (z) =ty(z) - ip v (z) , then R v (r) = r^ 1 (b^ + b +2 r + ...) 
l R v (r)|<. y(l + 1 r I / p+ |r 2 | /p 2 + ... ) jr^ 1 ) /p' Ha 

= Y I ^/P I V+1 / (1 - |r/p|) 

■ °( l r /p | V+1 ) . 

But c v /c v+i ’ r V r) 7 Vi (r) 

- r [ 1-R v (r) /Hr) ] / [l-R^Or) /^( r )j 
= r[ 1 + o (| r/p | V+ ' L ) ] 

and choosing an appropriate p, then | r/p f <6. In conclusion, since 
0 < 6 < 1 we may conclude that lim c^/c^-*- r * Also jr/p|is a common 
ratio of (c^/c^.^ - r) thus making the sequence geometric. (This fact is 
very essential before the application of Aitken's <5 process. ; which we 
will discuss later in this chapter.) This theorem may be extended 
to include the case in which the power series h(z) has p simple 
poles in [ z | <R , though in the examples that we have considered, this 
extension will be impraticable to apply. 

THEOREM B [5] Let h(z) be a meromorphic function in [z|<R having 
exactlv d Doles in - _ r not* npres.sari 1 v H-f Qf-f-npi- -fn t-Vip 

_L J £ 9 * ~ p J ~ — 


disk* Let 
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It-jj ^ tr^ $ ir^l ^ 4 •tolp < ° r R<’R and let 

•P(z) = 1 + ajz + ... + a p z p E (1 - r 1 1 z) ... (1 - r~ p z) 
setting K l ,(z) = K f ^(z)il +«(z) + ... + <rz p where 
k v p(z) is defined as H vp K„ p (z). Then 
= a i + o(v v ), and K v (z) = <p(z) + o (<r v ) 

this section, we will extend the coverage of Theorem A by including 
the case when p =» 2 . Let r^ and r 2 be 2 distinct roots. Then 

z) = V 5 (z)h(z) = b D + b lZ + b 2 z 2 + ... is analytic in |z| < R and 
we have 1 + a^z + a 2 z 2 +■...) (c Q + c-^z + c 2 z 2 +...) = 

A 

^ (z) = h Q + b^z + b 2 z + . . . and equating simular coefficients 

c = b 
o o 

C 1 + a l c o “ b l 

c 2 + a l c l + a 2 c o = b 2 


c v + a^c v + a2C v _2 = b v 

Suppose r represents r 1 or r 2< Multiplying the equations above 
by 1, r, r2, . ,.r v and adding, we obtain 

(1) r l c v + ( 1 + a^r)r v ^ c v _i “ ^ v (t). In a similar way we obtain 

(2) r v+1 c y+1 + ( 1 + a 1 r)rc v _ 1 =^ +1 (r) and 

( 3 ) r v+2 c v+2 + (1 + a i*')* V+ \.+2 = ^\ +2 (r) * 

and multiplying the first equation (1) by r 2 and (2) by r , the equations 
then show that the determinant 
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where they express the linear dependence of the three columns. 

Since ^v(r) = ^ (r) - R„(r), it can be seen that in the limit 

both r x and r 2 satisfy the equation K v2 (z) = 0 or equivalently, 
k v2 < z > = °* However , in case ^ = r 2 is a pole of h(z) of 
multiplicity 2, a further argument is necessary to show that the 
limiting equation has r^ as a double root. Accordingly, a 
different argument has been developed by Householder [5] in which 
confluence does not require special consideration. 

At this point, it is appropriate to develope the difference 
equation enabling us to compute the coefficients of the power series 
which was derived from the quotient g(z)/f (z) , f(z) *). Consider 
f(z) a Q + + ... which has a root r in its region of analyticity 

and g(z) fc> 0 + k^z + ... which is analytic within the same disk 
as f(z) but not vanishing at z = r, then 


g(z)/f (z) = h Q + + ... 

Therefore, 

g(z) - (a Q + a x z + a 2 z 2 + . . . ) (h 0 + l^z 2 + ...) 

= b o + b l z + b 2 z2 + • • • 

On equating similar coefficients, we have 



where ^ 0 


a n h = b 
o o o 


a h. + a,h * b, 
o 1 1 o 1 


a o h 2 + a l h l + a 2 h o * b 2 


a„hu + a.h, , + a 0 h „+...+ a h = b, 
o v i *v_i 2 v-2 n v-n 


4 

h 


making it possible to compute while making some; arbitrary se- 
lection for g where 1^ /h v+ ^ approaches r for some sufficiently > 
large V . 

. I 

With the knowledge gained from previous results we can \ 

compute a very useful expansion for approximating the root of least 

i 

modulus for a particular transcendental equation. This is known 
as the WHITXAKER-BERNOULLI expansion. In reference to theorem A, 
the root of least modulus may be expressed as 


r m h c /h^ + (hjL /^*2 ” ^o^ b l^ ^ b 2^ b 3 — b l^ b 2^ * * * 

= h Q /h 1 + (h^ - h 0 h 2 ^/h 1 h 2 + (h^ - h.^)/ h 2 h 3 + 

and futhermore, suppose h(z) - l/f(^) or h(z)f(z) = 1 where 

h(z) is meromorphic for | z |<R having a simple pole at z 3 r 

2 

within the disk. Let f(z)h(z) = h Q + h^z + h 2 z + ... 
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then multiplying and equating coefficients 
a o^o = ^ 

a o h l + a i h o = 0 • 

a o h 2 + a l h l + a 2 h o = 0 


a o h v + a l h v-l + ••• +a v h 0 = 0 

Solving the system of v +1 equations in v+ 1 unknowns we find that 
h o = 1 / a o 

h X = -aj/aJ 

h 2 “ a l t a l ~ a 2 /a o 


Thus we find r to be 




11 


Ocassionally , the sequence relulting from the application 

of the Bernoulli method to a particular transcendental equation may 

not provide satisfactory convergence to the root of least modulus, i.e. 

subsequent recursive computations of determines a sequence 

of slow convergence. To accelerate convergence, there exist a method that is 

2 

both highly effective and convenient, namely, Aitken s 6 method. 

A sufficient condition for the employment of this method is that the 
given sequence be geometric. Recall that this condition was established 
in Theorem A. 

Consider a sequence , a 2 , ... , ... converging geometri- 

cally to a limit, say, a, for some v , where a = h /h . Now 

v v v+i 

2 

the Aitken <$ algorithm may be derived in the following manner. Under 

the assumption that the sequence above converges geometrically to a 

V V 

value a or a-a » bh +e (h ) where 0 whenever h -HD which 
v v 

may be expressed as a - a^ = bh + o(h). From this result, we obtain 
(a-a v+1 )/(a-a v > = h + o(h) and (a -a v ) / (a-a^) = h + o(h) 
thus upon subtracting we have 

(a - a^ fl )/(a-a v ) - (a-a^) /(a - a^) = o(h) where 

2 2 V 
ViV - % - a( Vi ' + a ^i 5 ■ o(h 5 

Since 

a t - 2a + a -bh (1-h) 2 - o(h v ) thus 
v~l v v+1 

a = ( a v _! a %H . 1 - a v )/(a -i _2 V Vl ) + ° (hV+l) 



12 


In the following discussion, we will investigate a method 

i 

of solving transcendental equations whose theoretical developments 

are unrelated to the previously discussed materials. An algorithm 

will be developed which will enable us to locate the real root of 

least modulus through summing the reciprocated roots of a power 

series. The motivating reason underlying such an undertaking is 

that situations do exist where a transcendental equation may have 

a power series representation. 

Using basic principles from the theory of equations one can 

state explicitly the sum of the roots of a regular polynomial equation. 

n n“l 

We focus our attention to the general case where a Q x + a n _^x +. ..a o =0. 

Upon dividing through by a Q ^ 0 we have (a n /a 0 )x n + (a n _jya o )x n + . . .+1 =0 

2 n 

then 1 + CTjX + c^x + ... a n X = 0, In addition from the theory of equations, 

we know that the sum of the product of the reciprocals of ^ taken k at 

a time is (-l)a . It also follows that if given the equation 

k 

9 9n 

1 + a^x^ + ...+a n x =0 then the sum of the squares of the products 

k 

of the reciprocals of its positive roots taken k at a time is (- 1 ) (a^)* 
Suppose if given the real roots of a particular equation, say 

» ^ 2 * ^ 3 * where we are find the sum of the reciprocals of 

various powers of the roots. First, we may represent the sum of the 
reciprocals of the given roots as 

(1) 1/\ 1 + 1/X 2 + ... + 1/X n = a 

(2) 1/X 1 X 2 + 1/A 2 X 3 + ... + l/X n _ 1 X n - a 2 



Now assuming that Xi‘X 2 ’ A 3 * * * * are roots t0 the e< l uation 

ty ' 

1 + a ^x + a2 xZ + ••• = 0 and combining equations ( 1 ) and ( 2 ) we have 

(3) l/A* + H\\ + 1 A| + ... + 1/ A* = _ 2 « 2 

In a similar way we can find 

(4) 1/^ +1A 2 +IA 3 + + = 3/2(lA J + ... + 1A*) x 

(lA + lA 0 + • • • + i/a ) 

1 2 n 

+3(1 A iA ^ 3 + 1 A jA ^ ^ + • • 

-1/2 ( 1/A^ + I/A 2 + ••• + 1/A 

Thus from equations (1), (2) and (3) we have 

(5) 1 /A^ + 1/A 2 + ... + 1/A^ - 3 a i a 2 -3 a 3 and similarly 

n 4 

(6) £ 1 /A. = _ 4a,et„ + 2 a 2 + 4a a - 4a, 

v ' 1 k 1 12 2 1 3 4 

At this point it should be noted that the reason for seeking higher 
powers of 1/A ( k= 1,2,3,...) is to obtain a better approximation to 

K- 

the root of least modulus. Also in the manner in which (5) and ( 6 ) 
were derived, one may extend this method to include higher powers of 
1/A^ if further "refinement" is desired. Note that the rapidity of 
convergence of the expressions on the left of (5) and ( 6 ) is such that 
all but a finite number of terms will cluster about the "dominating" 
term thus allowing us to approximate the root of least modulus. 



CHAPTER III 


EXAMPLES OF SOLVING TRANSCENDENTAL EQUATIONS 

The methods of solving transcendental equations that were indi 
cated In Chapter II will now be applied to particular equations. The 
motivating reasons for the employment of these methods are consistent 
with those outlined in the abstract. 

Among the facilities used for aiding computations was an APL 
360 computer terminal. Because of the limited capacity of our compiler- 
some difficulties were encountered in cases where the coefficients of 
series under consideration converge absolutely at a fast rate. 

However, inspite of these limitations, considerable progress was made 
in making the objectives of this paper a reality. 



SUMMATION OF THE RECIPROCALS OF ROOTS 

Consider some transcendental equations and their solutions: 

(1) sin x /x = 0, using the Maclaurin series expansion for 

sin x, we find sin x / x = 1 - x^/6 + x^/ 120 - x^ / 5040 + ... 

From previous considerations, we see that = -1/6. 

a„ = 1/120, a„ = -1/5040. Now if we let P be a 
2 3 n 

root of sin x/x = 0, then the sum of the reciprocated 
fourth powers of the roots would be represented by 
Il/^[=l/90, thus Pfj/— /90«3.14. The actual value is 3.14+ . 

Of course, an even more accurate approximation may be 
made by calculating the sum of the reciprocated eighth 
powers of the roots. 

(2) Consider J Q (x) = 0 where J Q (x) is a particular solution to the 
Bessel equation. J Q (x) = 1 - x^/2^ + x^/2^4^ - x^/2^4^6^ + ... = 

ct = - 1/4, a 2 = 1/64 , 0 I 3 = - 1/2304 

Letting P^ be the n^ positive root of J Q (x), thus we have 

21/P^ = 1/4 where El/P* = 1/32 and El/p^ = 1/192 

thus P = *Vl92 ~ 2.40. The actual value of this root is 2.4048. 
n 

(3) In our investigation of tan x » x, and tan x, and tan x being an 
non-integral function, we may write an equivalent expression 
for the given. Namely, sin x - x cos x = 0, cos x / 0. Using 
the power series for cos x and sin x 
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we have the following series: 

x^/3 - x^/30 + x^/840 - x^/45360 + . . . = 0. Multiplying 

3 

both sides of the expression by 3/x , we have 
1 - x^/10 + x 4 / 280 - x^/15120 + ... = 0 ; x ^ 0 . 

We see that = -1/10, <* 2 = 1/280, « 3 = -1/15120 

th 

assuming that P n is the n positive root of 
sin x - x cos x = 0 and from previous results , we have 
Sl/P^ = -1/10; Sl/p 4 = 1/360, 21/p|j = 1/7875 

g 

thus 75 ^ 4.46. The actual value of this root is 4.4934, 

(4) Now consider another equation which is associated with the 

trans-vibration of beams, namely cosh x cos x =^1 . Again utilizing 

the Maclaurin series expansion for cos x and cosh x the expression 

becomes 1 - x 4 /12 + x®/5040 - ... =0 

2 1/P 4 = 1/12 and P = \/ 12 =1.86 
n n 

also 2 1/P® = 11/1680 and ? n = yj 1680/11 = 1.87 + ♦ The actual value is 1.8751. 

THE METHOD OF WHITTAKER - BERNOULLI 

We now focus our attention on another method which may be 
alternately used to approximate the roots of a transcendental equation. 

This method is known mainly by the Whittaker-Bernoulli method which 
was derived in chapter II. For some equation say 

x + » » . =0 the smallest root can be approximated in 

absolute value by the given expression: 
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a 2 a 3 

a l a 2 


X = a o /a l “ a o a 2 /a l 

a l a 2 

- a 

0 

a l a 2 


a l a 2 a 3 


a o a l 


a o a l 


a o a l a 2 


1° a o a i| 

Now we employ this method to treat the equations introduced previously. 

= 1 - x /2 2 + x^/2^4 2 - x^/2^4 2 6 2 + ... where we let z = x 2 
then J c (x) = 1 - z/2 2 + z 2 /2 2 4 2 - z 3 /2 2 4 2 6 2 + ... - 0 
and 


1 


3 


1/64 1/2304 


z = 1/1/4 + 


l 2 1/64 

-1/4 

1/64 

-1/4 1/64 


“ 

-1/4 

1/64 

1 1/4 



1 

-1/4 


-1/4 1/64 1/2304 

1 -1/4 1/64 

0 1 -1/4 


z=4 + .3333 + .16447 


thus x = 2.40 


Evaluating the equation cosh x cos x = -1 

Zj, ry 

where x = z and cos z cosh z = 1 - z/12 + z^/5040 - ... = 0 
and a- = 1, a x = -1/2, a £ = 1/5040 


then the root 

z = 


1 /i / r rs f /\\ 

. L * 


1.8744 


-1/12 ,1/5040 

1 . 1/12 
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Consider the equation sin x/x » 0 whereby applying the Maclaurin 

ii 

series we have sin x/x = 1 - x^/6 + x^/120 - x^/5040 + ... =0 

where a Q ■ 1, a - -1/6 , a 2 = 1/120, a 3 = -1/5040 and letting 
2 

z- x 


so z = -1/1/6 + 


1/120 -1/5040 


1 - 1/20 

■ - 1/6 1/120 

1 


- 1/6 1/120 
6 j 


- 1/6 1/120 


- 1/6 

1/120 

- 1 / 504 C 

1 - 1/6 


1 - 1/6 


1 

- 1/6 

1/120 



0 

1 

1/6 

+ .9207885 

or x = 

V 9 . 49279 

= 3 . 

081 



Now let us consider the equation 

tan x = x or in power series form 
1 - x^/10 + x^/280 - x^/15120 + . =0 and letting z = x^ 

we have 1 - z/10 + z^/280 - z^/15120 + • . . = 0 
where a = 1 , a. = -1/10 , a = 1/280 , a = -1/15120 

o A. 2 J 


and x - 4.49 



THE EXTENSION OF 
BERNOULLI’S METHOD 

l 

Consider the Bessel function J Q (x). Following closely the 
theorem outlined in chapter 2, we shall obtain an approximation to 
the root of : 

J D (x) = 1 - x 2 /2 2 + x^*/2 2 4 2 - x^/ 2 2 4 2 6 2 + ... ■ 0, and letting z =x 2 
J c (z) = 1 - z/2 2 + z 2 /2 2 4 2 - z 3 /2 2 4 2 6 2 + ... = 0 

we find that h Q /h^ = 1, h-^/tv? = 5.333, h 2 /h 3 - 5.684 , h 3 /h 4 = 5.763 
h 4 /h 5 = 5.779, h 5 /h 6 = 5.782, h 6 /h ? - 5.783 

Noting the relative magnitudes of the last 2 terms in the sequence, it 
becomes apparent that the sequence has begun to converge to a root. 

Thus we find that z = 5,783 or x = 2.405 

Now we consider a function which converges somewhat less slowly, namely 
log(l + z) = 1/2. Expanding the function by the Maclaurin series, we 
have -1 + 2z - z 2 + 2/3 z 3 - 1/2 z^ + . . . = 0 

Again computing recursively, h/s from the difference equation 
h v = 2 ( h^ + l/2h y _ 2 + 1/3 h + ...) and h Q =l, 1^=2, h 2 =3, h 3 =14/3 
h 4 = 43/6, h 5 = 166/15 , h 6 = 767/45 and 1^/h = .5, h 1 /h 2 = .667 
h 2 /h 3 » .643, h 3 /h 4 = .651, h 4 /h 5 = .648, h 5 /hg = .6493. According 
to this method, the indicated root is .6493. The actual root to the 
above equation is e"^ 3 - 1 = .64872. At this point it would be inter- 
esting to see whether the Aitken 5 process [5] would have given us 
better accuracy while using fewer terms. This process is stated 

aS X ~ X n+2 " (x n+2 ' x n +l )2 7 (x n +2 ~ 2x n+l + V ’ thuS 

we find that x » .6493 - (.6493 - .648) 2 /(.6493 - 1.298 + .651)=. 6486. 

The result compares very favorably with the actual least root of 
the above equation, namely , log(l + z) = 1/2. 



Let us now focus our attention on the equation - tanv4/v£ « 

-1 - 1/3 z - 2/15 z 2 - 17/315 z 3 - 62/2835 z 4 - 1382/155925 z 5 - 

21844/6081075 z 6 = 0 and computing h Q , , h 2 , .... , h^, which are 

respectively 1, -1/3, -1/45, -2/945, -1/4725, -2/93555, -1382/648750375. 

The ratios h Q /h^ , h^/l^, ..., h^/h^ becomes respectively 

-3, 15, 10, 5, 10 , 9.9, 9.870. On the other hand, from the 
2 

Aitkin 8 process we find that the indicated root is 9.938. It 
should be noted here that we used only four (4) terms from the sequence 
to make this evaluation! 



GRAEFFE' S METHOD 


In computing the root* of least modulus to the equation 
log(l + x) =1/2 whose power series representation is 
-1.+ 2x - x 2 + 2/3x 3 - l/2x 4 + 2/5x 5 - l/3x 6 + 2/7x ? .... we find 
that the root of least modulus is .6448 where the actual root to 
this equation is .64872. Eight (8) terms and iterations were used. 

For equation (2), we have sin x/x = 0 where x f 0 and 
sin x/x = 1 - l/6x 2 + l/120x 4 - l/5040x^ ... Employing the Graeffe 
Method to these coefficients, it was determined that root of least 
modulus x = 3.078 + ... On the other hand, the actual value of such 
root should have been 3.14 + ... Four (4) terms and iterations were 
used . 

In applying the Graeffe Method to the Bessel equation of order 
zero, namely J o (x) = 0, to the coefficients of the Macluarin series 
(which was stated early) we find that the root of least modulus is 
x = 2.567 + where the actual root is x = 2.40 +. • Six (6) terms and 
6 iterations were used. 

In reference to the equation cos x cosh x = -1, it was deter- 
mined from the application of the method of Graeffe, that the root of 
least modulus is x = 1.875 where the actual value of this root is 
given as 1.87 +. Three (3) terms and 6 iterations were used. 

Finally, in computing the roots to tan x/x = 0 and with 6 
applications of the Graeffe Process, we find that the root as indi- 
cated by computation to be x = 4.147 where the actual root is 
x = 4.49 +. Four (4) terms and 6 iterations were used. 

See Appendix II for the actual computation of the roots. 



CHAPTER IV 


REMARKS AND CONCLUSIONS 

In Chapter III, we mentioned various existing methods of 
globally approximating the roots of transcendental equations. The 
various approaches mentioned are advantageous from a particular stand- 
point depending upon the type of solution desired, size of the root 
sought, computational aids available and more importantly, the effect 
of "roundoff" effor on the approximation. At this point, we discuss 
the suitability of the various methods explored in this paper. 

The method whch employs the summing of the reciprocated powers 
of the roots of the transcendental equation offers fast convergence. 

Implicit in this method of approach is that one may obtain a 
fairly high degree of accuracy with the least number of terms 
from the power series representing the equation. In economizing 
the terms of the power series we reduce the effect of roundoff error 
on the root approximation. In addition, this method allows one to 
achieve good accuracy, especially, if the computational aids are 
limited only to the slide rule or desk calculator. 

On the other hand, the Whittaker-Bemoulli method offers a 
somewhat different approach as an alternate to approximating the 
roots of transcendental equations. This method allows one to 
obtain a relationship between the root of the transcendental equation 
and the coefficients of the truncated power series. One advantage 
here is that we can reduce C.P.U. time by as much as 20% over the 
Graeffe method, and at the same time, this method yields acceptable accuracy 



(for most cases) with the least amount of effort. Again we reduce the 
effect of "roundoff" error by employing fewer number of terms of the 
series. 

The method of Bernoulli [Koenig's extension] does guarantee 
convergence to a root but not nearly as fast as other methods as 
indicated in the application to J Q (x) = 0. One of the advantages of 
this slower convergence is that it may provide suitable separation 
when the roots are nearly of the same magnitude. However, in the 
case of J Q (x) - 0, convergence was accelerated by employing the 
Aitken 8 process. This was possible because h,./h r ^ ( 0,1,...) 

converges geometrically as was shown in Chapter II. In using the 
Bernoulli method alone, it sometimes requires a large v inorder to 
obtain accuracy comparable . to the other methods mentioned in this 
paper. 

Among the methods that have been studied in this paper, the 
Graeffe method is the most "global" in that it is capable of yielding 
all the roots to a transcendental equation from a theoretical 
standpoint. On the other hand, there are some inherent drawbacks 
associated with the process when applied to transcendental 
equations. First, when the power series is truncated and the Graeffe 
method applied, some sign alternation may result. When this 
situation is observed in the case of polynomials, then one is 
immediately alerted of the presence of complex roots. On the contrary, 
truncating the power series representing a particular transcendental 
equation tends to obscure the properties or behavior of roots 
subsequent to the first. Throughout the examples that were 



1,lcluded (aee rI) * «■> —ti. behavior was not evldent 

when the root of least modulus was computed from the first coefficient 

I 6 r d - - roots subsets to 

rSt “ he ” COnPUted fr ° n thC re “i-in S coefficients. Secondly 

situations were encountered where the coeffiM. , 

• . , „ coefficients would "absolutely" 

o zero for increasing subscripts. In attemnf 

* attempting to compute 

more terms to refine rto „ 

the accuracy in determining other roots and 

iterating Craeffe's process, one is rewarded with an inordinate 

number of zero coefficients m.hing it nearly possible r„ determine 

nformation about the roots. This is especialiy evident in the 

example displayed on page 31 of the appendix. In a positive 

sense, in observing rbe examples of Chapter Hi, it can be seen 

that the method of Graeffe yields acceptable accuracy and fast 

convergence in computing the root of ieast modulus. Also this 

method seemed very comparabie with the capabiiity of the API 

programing ianguuge. Computer programs for the methods mentioned 
in this paper are available in Appendix I. 

In the case involving polynomials. Bernoulli's method may 

‘ ‘ PPlIed Pri ° r C ° -*-* Inorder to increase the sep- 

aration of roots [1J. Since the domain of our investigation did 

not include the approximation of real roots nf i 

real roots of larger modulus, it 

-P be interesting to h„ow„ the degree of separation of toots 
that can be expected from the employment of the above procedure 
to a truncated power series. As a contribution to this effort. 

to what extent shouid we incorporate the method of asymptotic 

Expansions to the power serine r* 

P Serl ° S presenting a particular trans- 
cendental equation? 
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APPENDIX I 


APL COMPUTER PROGRAMS 



27 


A 
A l 


PL PROGRAMS RELATING TO TRANSCENDENTAL EQUATIONS 


7£>ET[p]V (WHITTAKER -BERNOULLI METHOD) 

V Z+A DET B ',1 ;J ;A' 

Cl] 1 O +/1 ) p ( , ^ ( A ,.4 ) p3 ) , \ I*A. , K* \ 0 

[2] ^(V/0=M[1;1) J)pO 

C3] -K 0 *l<-r-pppj?«el<}>( j -,1 ):)>[1] B-(Jx3t ;1] )*,xZ?n ;]-*-.U[l ; ] 1] )/2 

[4] Z+x/Z 
7 


ERECT R] V (METHOD OF BERNOULLI VIA KOENIG’S EXTENSION) 

V Z+RVC B 
Cl] S^-ppFt-,<|>5 

[2] '+2x(pB)>pZ*-Z ,P + ,x(-pB)iZ 

V 


N COLUMN (GRAEFFE’S METHOD) 

NU1+COE V 

[1] COUNT * 0 

[2] 01«-(l.(7?-l)p ~2 2) + .*(((O t -2\pV)*U) AU<rT° .<)*(V*2) "l o.x"l + r^ u p) 

<t>(2,(~l+;?^ro,5xp7),pF)pF, T / 

[3] COUNT*(COVNT + 1) 

[4] O2*(70tf5!r 

[5] HCOUNT -K )/ 0 

[6] 

C 7 ] -*-2 

[8] N+l 7 



APPENDIX II 


EQUATIONS EVALUATED BY GRAEFFE’S METHOD 



u 


COE SH ^ cos x CwSft X = 0) 

2 

1 0.006547619048 3. 936759889ff“8 

4 

1 4.279257999F“5 1 . 5 49807842E“l 5 

Z 

t 1 .8 31 2 01 8 032? ""9 2.4019O4349£“3O 

16 

i 3.3533OOO40ff“i8 5.769144SF~60 

i2 

1 1 , 1244621 1&F~3 5 0 

64 

1 1.264415055F”70 0 

128 

10 0 
256 

1 ' 0 0 


R+* 1. 26441 5055F“70*< *256) 

R 

1.875165023 



COE XS (Tan X /X ■ 0) 

2 

1 0.002857142857 ■'4.724111867£~7 4 

4 

1 9 . 108 0876 B2?~6 “2 . 47 721 285 5i?“ll 1 

8 

1 1.32501 51 B32f*’l0 2 . 6 512 04 342/T 22 

15 

1 1 ,702641148J?'*20 ~2.672562663£“44 

32 

1 2.899S21391F~40 0 0 

64 

10 0 0 
128 

10 0 0 
256 

10 0 0 


i?«-U . 702641146£’20*< *32 ) 

R 

4.147414124 


.374177655ff"9 
.91334301 5J?~17 
3.6608 8 1 49 5£~34 
1.340205332£"67 



& 


COE JX (J o U0 = 0) 

2 

1 0.03125 0.0001899401347 1. 094S972l9E“8 6 .768187595E"16 1 . 321996248£“26 

4 

1 0, 0005966822306 3. 5 39313286ff"8 1 . 195S719722?“l6 4 . 577942226ff“31 1 . 747674078ff"S2 
8 

1 2. 8 52 4341 8 65** 7 1 . 2 52531 178£~15 1 .4293890992?“32 2 . 09575 5085F”61 0 
16 

1 7.B85874552F"14 1 . 5 68834 344J?”30 2 . 0431 5 31 96£~64 0 0 

32 

1 6.215564077JF"27 2.46l241198j?“6 0 0 0 0 

64 

1 3.863323187J?“S3 0 0 0 0 

128 

1 0 0 0 0 0 

256 

1 0 0 0 0 0 


J?«-*(6 . 21 5564077F~27 )*( *64 ) 

R 

2.567301452 



c»> 


(70F SX ^SIN X /X = Oj 

2 

1 0.01111111111 3. 306878307£~6 3 . 936759889ff“8 

4 

1 0.0001168430335 “8 . 639000868£“l0 1 ,549807842K“l5 

8 

1 l.S38009465ff“8 3. 841 548606F“19 2.4 Ci904349E' 30 

16 

1 2 . 35779 0018F~16 7. 3%9192446£'”38 5 ,76§1445£”60 

32 

1 5. 5591590 32 F*" 32 2. 7l 001 3468£~75 0 

64 

1 3.090424914F"63 0 t> 

128 

10 0 0 
256 

10 0 0 
P 


J?«-*3.09Q4249i4g"63*C*128) 

R 

3.078642304 



